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1 Introduction 



New RHIC data on the multiplicity in Au+Au collisions |J and its centrality de- 
pendence P, [H has given us new insight into the dynamics of ultrarelativistic heavy 
ion collisions. The data from central and nearly central collisions can be understood 
@, H in terms of a conventional soft + hard two-component picture |J, but also in a 
dynamically more unified perturbative QCD + saturation model [7], ||. 

The results of @ are formulated in the form of scaling rules: quantity ~ CA a (yfs) b 1 
where the constants C, a, b are determined numerically for central A+A collisions using 



independently determined parton distribution functions || with shadowing fL0 |. For 
example, for the dominant saturation scale p sat and for the multiplicity A" per unit 
rapidity one finds: 

p sat = 0.208 A ' 128 v^ ' 191 , (1) 
N(p sat ) = 1.383 A ' 922 v^°- 383 , (2) 

where p sat and \/s are in units of GeV. Thus the multiplicity A" ~ ^40.922 j ns t eac i Q f 
the exponent A" ~ A 4//3 appropriate for hard collisions or A" ~ A appropriate for the 
saturation model if xg(x) ~ const and a s ~ const. Even more striking is the rather 
fast powerlike dependence A" ~ i/i ' 383 , much faster than the ~ log(s) or ~ log 2 (s) 
behaviour observed for pp collisions. 

The results in ([[])- (0) lead to rather definite and easily testable predictions for the 
overall magnitude and the A and s/s dependencies of heavy ion experimental results. 
Agreement with the first RHIC results for the charged multiplicity and -^i dependence 
is good || . It would thus be of some value to derive the numerically computed scaling 
parameters analytically and to understand the underlying physics. It is the purpose of 
this note to carry out this derivation. 

The dominant part of the analytic estimate is a derivation of an accurate approxima- 
tion for hard production of minijets in pp collisions. This, of course, is a most standard 
problem, but with one important difference: as we are interested in the physics of heavy 
ion collisions in the RHIC-LHC energy range, y/s > 100 GeV, and for A ~ 200, we know 
the magnitudes of the dominant scales of the problem, Q ~ p SSut and x ~ p sa t/y/s. It is 
thus sufficient to find out an approximation for the gluon distribution functions in this 
range. In fact, we shall find that the very simple estimate xg(x,Q 2 ) ~ (Q/x) s with 
5 ~ 0.5 is quite accurate in the vicinity of the dominant scale Q, which depends on x. 
After obtaining an analytic approximation for the perturbative minijet cross section 
da/dy(y = 0,Pt > Po) (Section 2) it is straightforward to find the scaling exponents 
(Section 3). 



The idea of saturation originates from [O, O, [T3| and is usually discussed in con- 



nection with small-x behaviour of gluon distribution functions ("initial state satura- 
tion"). Here we have in mind rather a picture with the saturation of produced gluons 
("final state saturation"). Using the analytic approximation we can show the close 
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phenomenological relation of initial and final state saturation (Section 4). This close 
relation also follows from the fact that N(p SSLt ) is proportional to the initial state gluon 
distribution probed at the final state saturation scale. The analytic approximation can 
also be extended to models of local saturation (Section 5). 



2 Analytic estimates of the minijet cross section 

Consider first inclusive gluon production from the subprocess gg— >gg in pp collisions: 

da „ r , , . , n,9a? / x 2 x " 



(l!/(Pi>j K J d V2 x ig {x^%) X2g ^ P %)^[l-^] . (3) 



Here K describes the effect of higher order corrections |14j] and the fractional momenta 
are X\ = \xT(e y + e y ' 2 ), x 2 = \xT[e~ y + e~ m ) with xt = 2pt/\/~s- The integral is 
over — log(2/xr — e~ y ) < y 2 < log(2/x r — e y ). The last factor in (Q) is equal to 
(1 + 2 cosh Yf /{2 + 2 cosh Y) 3 , Y = y - y 2 . 

For the analytic estimates, we need to approximate the gluon distribution xg(x, Q 2 ) 
in a region which dominates the Pt- and ^-integrations. From ([2J), we see that p sat ~ 1 
GeV at v/i = 200 GeV, and p sat ~ 2 GeV at ^/s = 5500 GeV for A ~ 200. In addition, 
as shown in Fig. 1 of J/]], we know that about 90% of all the minijets produced above 
the saturation scale have transverse momenta p sat < Pt 2p sa t- The scale in (^) 
is chosen as Q = pr, so the relevant region in Q thus is 1 ... 2 GeV for y/s = 200 
GeV and 2 ... 4 GeV for = 5500 GeV. In what follows, we shall fix y = in (§. 
This makes the integral even and it suffices to consider only the integration region 
< y 2 < log(2/xT — 1). The fractional momenta are now limited to xt/2 < X\ < 1 
and xt/2 ~ (2/xy — < x 2 < Xt/2. From the numerical computation of @, we 
have checked that the dominant (> 70%) contribution comes from y 2 < 4 for -^i = 5500 
GeV (p T > 2 GeV), and from y 2 < 3 for yfe = 200 GeV (p T > 1 GeV). In the analytic 
estimates we thus need to describe xg(x, Q 2 ) at Q/\/s <x< 50Q/y/s in the vicinity of 
Q = 2 GeV for yfe = 5500 GeV and at Q/^fs <x< 20Q/v^ in the vicinity of Q = 1 
GeV for ^ = 200 GeV. 

The gluon distribution xg(x, Q 2 ) as obtained from the set GRV94-LO in the x, Q 
region discussed above is shown in Fig. la. The symbols are for ^Js = 200 GeV (circles) 
and for y/s = 5500 GeV (squares). For our purposes, and as the kinematic region is 
now limited, a simple power law x~ 5 , shown by the solid lines in Fig. la, reproduces 
xg adequately. To simulate the effect of pQCD scale evolution, we note from Fig. lb 
that in the dominant region (to the right of the thick tilted lines), xg(x ~ Pt/\/s,Pt) 
is approximately constant. This suggests a simple fit 

xg(x,Q 2 )=C (^^f } (4) 

x 
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Figure 1: (a) The gluon distribution xg(x, Q 2 ) as a function of x at fixed scales Q, plotted in 
the region which dominates minijet production at saturation at -^/s = 200. ..5500 GeV, i.e. at 
x ttQ/ y/s with Q = 1, 2 and 4 GeV (see the text for details) . The squares and circles show 
the GRV94-LO distributions §§, the solid lines are the fit xg(x,Q 2 ) = 0.32(Q/GeV/x) a5 . 
The solid thick line shows the gluon density probed at saturation, xg{xT S a,t,P 2 & t) with A = 208 
at y/s = 200 . . . 5500 GeV. (b) The gluon distribution xg(xT,Pr) as a function of the scale px 
for \fs = 200, 500, 1500 and 5500 GeV. The symbols are the GRV distributions and the thin 
solid lines are the fit. The thick curves (solid, dotted, dashed) are the gluon densities probed 
at saturation, xg(xT S a,t,Psa,t)i a t different \fs for A = 208, 107 and 40, correspondingly. The 
region relevant for our problem is to the right of the thick curves. 



with 5 = 0.5 and Cq = 0.32, as shown by the solid lines in Fig.l. Below we shall see that 
iVfeat) ~ xg(x Ts ^p 2 J, with x Tsat = 2p sat /^ = 0AlQA 0128 (^/GeV)-°- 809 , where 
the relation ([I]) is used. The gluon density probed at final saturation, xg(xT sa ,t,Psa,t) * s 
shown by the thick line in Fig. la for A = 208 from y/s = 200 GeV to 5500 GeV. Note 
that on this curve each point in x now corresponds to one yfs and one p sat . 

If a similar procedure is carried out for the CTEQ5 set of parton distribution func- 
tions fl5|], effectively the same result is obtained, the value of Cq is decreased by less 
than 10% and S ~ 0.47 is somewhat smaller. 

With the fit xg(x,p 2 T ) = C (p T / GeV / x) 5 , (5 = 0.5, C =0.32), Eq. (§ can now be 
expressed as 

j „ ,/ 2\ I- 3 o i log(VI/PT-l) 
d(T o ^ 9 «s(Pt)^ / s 5l2 ^ n f3\A, n+ 6 f , . r i. 3/2, 



As our focus is at a region where Xt *C 1, we write 

10g(-v/s/PT-l) 

/ o ^ 2 [cosh|]- 2A = f dy 2 [cosh^]- 2A - f ^ 2 [cosh|]- 2A (6) 

log(Vs/PT-l) 
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1 rxx 7 

= B(A,i)-2 i / o dzz"- l (l-^ A -' (7) 

- ^.^(^grr)^-^ (« 

where B(A, |) is the beta function, A = n + 5, and in the second integral a change 
of integration variable from yi to z — xt/x\ has been made. In the limit of small xt, 
the leading term for each n is given by the first term with the beta function. The pt 
distribution thus becomes 

d ° ~ 2 ■ K ^f [C ( ^ /GeV )rB(5, l)a(J) + Q(4) 0) 



dyd 2 p T 2p% 1 uv 2 7 J v ' 2 

where 

<l) = 1 ** + * ^tl) JL ^Mff ^ 0-6904 (10) 
4± + 5 4 2 (± + £)(§ + 5) 4 3 (± + cT)(§ + £)(§ + 5) 

and 

1 ygXg g = o, 5 

2 } = f(rT5) = 77 

Notice that now, in (|), [C (^§^) s ] 2 = [xg(x T , p 2 T )} 2 appears. 
Integration over then gives the minijet cross section 



, . da 
2a pQC »(p ) ~ - 



- , 27rp T dp T -^- (12) 
»=° Jpo dyd 2 p T 

« 2-^E( ( 5,i)a(5)N(x r o,^)] 2 X(^), (13) 



where 



•'Po IPtJ Po a qcd 



log (log 



A2 A2 "n*" to A2 

iv QCD ii QCD ii QCD 



(14) 



where E\{z) is the exponential integral and the running coupling a s (Q 2 ) is that to one 
loop. Denoting z = log(po/AQ CD ) and using the approximation |1| 

ze ^(z) = + e(s) 15 

z l + M + o 2 

where \e{z)\ < 5-10~ 5 , a x = 2.334733, a 2 = 0.250621, b x = 3.330657 and b 2 = 1.681534, 
we arrive at the following expression for the hard cross section at central rapidity: 



da 
dy 



Q7r 1 n 2 ry 2 (r> 2 ) 

^K-B(5,-)a(5)-f(\og^)^[xg(x T0 ,p 2 )] 2 , (16) 

p Z Z A QCD Po 



where f(z) = (0.995924^ + 1.430913)/(z + 3.330657 + 1.681534/z), x T0 = 2p /^ and 
xg(xTo,Po) = Co( v ^ // ^ eV ) <5 . Setting f(z) = 1 would correspond to the approximation 
a s{Pr) = a (Po) m the integral in Eq. (|I~4T). 

As anticipated based on the precision of the rough fit to xg, our analytic estimate (JTEf) 
reproduces the "exact" numerical result (with gluons only, no shadowing) to about 10% 
accuracy near po = 2 GeV at y/s = 5500 GeV. An improved accuracy would require a 
better fit to xg(x, Q 2 ) in the regions of large y 2 (large X\). The terms 0{J xto/2), now 



neglected, contribute only at the level of a few percent at yfs = 5500 GeV, po = 2 GeV, 
and at a level of 10% at y/s = 200 GeV, po = 1 GeV. Since the main emphasis here is 
to understand the origin of the scaling exponents, we leave the overall normalization 
as a rough estimate. 

To extract the A and yfs scaling exponents for the saturation scale analytically, we 
introduce a second parameter £ by noting that the complicated po dependence of the 
product a 2 s f in flilf) can, in the relevant range, to good accuracy be represented by a 
power: 

a-(Po) v / /( 1 °g(Po/A 2 QCD)) « £>(A QC dM^, (17) 

where D = 0.775 and £ = 0.444 (with A QCD = 0.2 GeV and N f = 4). The accuracy of 
this approximation is within 1.5 % in the region p = 1 . . . 2 GeV. 

Also the first py-moment of the j>r-distribution can be computed by using the same 
sequence of approximations as above. The result is 



d<T f 00 da 



« 2K^B(S,l)a(8).2f(log-^-)^-[xg(x T0 ,pl)} 2 , (19) 
^ ^ Aqcd Po 

where now the same function / appears as in Eq. (|1~6|) but with a different argument. 
For the average p?, we thus get 



(Pt) 



2/(log(po/A QC D)) _ „ ( Po 



y=o /(log(^/A^ CD )) V Aqcd 



■ftwF- r^ "Po (20) 



where in the last step the power law approximation again holds in the region p — 1 . . .2 
GeV and rj = 0.0624 and F = 1.399. 



3 The scaling exponents 



We can now apply the analytic approximations ([TBI) and (|17D to the minijet cross 

for central A+A collisions: 



section (|lj) in the final state saturation condition 

N(p ) = T AA (b = 0)2a pQC B(po) 



pIRI 



(21) 
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Saturation is a dynamic phenomenon and, in the weak coupling limit, there would be 
powers of a s together with various numerical constants in ([H]). Taking a constant value 
~ 0.3 for a s the net effect in is an overall constant of about 1. Even at the LHC 
one is most likely far from the weak coupling region and we shall not keep the coupling 
constant dependence in the right hand side of (|2l"l) explicitly. This approximation is in 
agreement with RHIC data. Note, however, that a s is kept in flUf). 

Using T AA (0) = A 2 /(7rR 2 A ) with R A « 1.12A 1 / 3 and Eq. ([TBD in the power-law 
approximation (0), one finds that the solution of ( |2~1"D is 



Psat 



3y/ B(6, i)a(6) 1 



•C n (~) s -D-Ai 



1.12 2 ■ (fm-GeV) 2 °V QCD 



V(2+0 



^1/(4+20^1/(6+30 ^/P+Q (22 



where now the origin of each factor can be easily traced down. The exponent 5 comes 
from the behaviour of the gluon structure function in Eq. (Q) whereas the exponent £ 
originates from the running of the strong coupling constant in Eq. (|17|). The numerical 
value of the constant in front of the i^-factor is 0.1625. It is also understood that p sat , 
Aqcd and y/s are in units of GeV. We have also kept K separate to show how p sat 
and, especially, iV(p sat ) depend on it. The initial multiplicity of produced gluons at 
saturation, iV(p sat ) = Pg at R A , then is 

N(p sat ) « 0.850 ■ 1^/(2+0^(6+20/(6+30 ^/(2+0 (23) 

Note that the dependence on the K factor is rather weak, N(p sat ) ~ K 0A1 - instead 
of ~ K. Substituting the numerical values for the coefficients Co and D, and for the 
exponents 5 and £ as discussed above, and K = 2 as in 0, we obtain the following 
scaling laws: 

p sat » 0.187A°- 136 v^ a205 , (24) 
N(p sat ) « 1.13A - 939 v^ a409 . (25) 

Since the numerical results in Eqs. ([!]) and (fj) contain shadowing, which is not 
included in the analytic estimates above, we should compare the scaling laws obtained 
above with the ones obtained numerically without shadowing (all parton flavours in- 
cluded): 

p sat = 0.193 A - 137 y/s°' 204 , (26) 
N(p sat ) = 1.20 A ' 941 a/s ' 408 . (27) 

The agreement is good and we have thus analytically understood how these scaling 
laws arise. 

The numerical result for the CTEQ5 set |15| (no shadowing) is p sat = 0.208 A 0A4:1 y/s°' 192 , 



-^(Psat) — 1.37 A°' 952 y / s 0,384 . The somewhat slower yfs dependence follows from a 



somewhat slower evolution in this set, S ss 0.47. 
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Based on Eqs. fl2~l|) and (|T6| ) the multiplicity of produced gluons at saturation can 
also be cast in the form 



N( Psat ) « ^K9B(5,i)a(5)^f(\og(plJAl CD ))aMJAxg(x Ts ^ P L t ) (28) 

Thus we see that the initial multiplicity of produced gluons directly probes the gluon 
distribution at the saturation scale, as derived in for initial state saturation. The 
powers of a s are not the same because they differ already in the saturation condition 

Nuclear shadowing effects can also be discussed in the analytic approximation. Over- 
all they are a fairly small correction to the results above: the numerical evaluation of 



N(p sa ,t) with the EKS98 shadowing pi shows a 16 % reduction at ^fs = 5500 GeV 



and a 7 % reduction at ^fs = 200 GeV for A = 208. For smaller nuclei the effects 
are smaller. Shadowing obviously slightly decreases the effective exponent 5 in an A- 



dependent way. The dependence of the factor J B(5 } ±)a(5) on A remains, however 



small. Disregarding the few percent effects from the factor ^/(log (p 2 at /AQ CD )), we 
arrive at the following simple scaling for the multiplicity of produced gluons at satu- 
ration 

iVfeat) ~ Aa s (p 2 aat )xg A (x Tsat ,p 2 aa:t ), (29) 

where now p sat is from Eq. (|l|) and xg A is the shadowed gluon distribution per nucleon. 
This result is tested against a full calculation of Eq. (|2|) in Fig. 3. The agreement with 
the numerically obtained results is good in the scalings with both yfs and A, especially 
at large \fs and large A(] If the initial state multiplicity is directly proportional to the 
final state multiplicity, the measured charged particle multiplicity then directly probes 
the nuclear gluon distributions at the (final state) saturation scale. 

Using the power- law approximation of Eq. (0), we obtain the average initial trans- 
verse energy per produced particle 

£r(Psat) T aa (0)<j pQ cd(E t ) 



(Pt) 



V=0,Psa,t 



1.546(|^)" - Psa , (30) 



-W(Psat) T AA (0)2a p QCD 
From here, using the analytic approximation for p sat from Eq. (|22|) and N(p sat ) = 
P 2 a t R % w e get 

E T (p s ,t) = 0.191ir( 3+ ")/( 4+2 ^A 2 / 3+ ( 3+ ")/( 6+3 «)v^ 5(3+,?)/(2+f) (31) 
= O.igiK - 627 ^ 1 - 08 ^ - 626 . (32) 
At saturation, the initial number and energy densities become 



iVGPsat) _ 1 3 _ E T ( Psat ) _ 4 1.546 ( Psat \ 
_ _ , , — X),, „ ^ « €.i — , , — T)„ I _ _ I 



0.0624 



n * ~~ 777 7 ~~ Psaf e « — T77 T ~~ Ps&t 7TT7 > v KPsax) — nn A /p : 



V(p Bat ) vr™ 1 V(p sat ) " sat vr VGeV; 



sat- 



(33) 



4 The slight kink in the curves with shadowing originates from taking the shadowing to be scale 
independent at p sa t smaller than the minimum Q in the EKS98 parametrization. 
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Figure 2: The initial multiplicity of produced gluons N(p SSu t) at saturation as a function of 
y/s. The symbols are the numerically obtained fits, Eq. (Q) with shadowing included (filled 
symbols), and Eq. ( |27| ) with no shadowing (open symbols). Circles, boxes and diamonds 
stand for A = 208, 107 and 40, correspondingly. The solid curves are the prediction of Eq. 
(|29|), withp sat computed from the numerical fits, Eqs. ([[]) and (|26|) and with GRV94LO gluon 
densities with and without shadowing. The solid curves are normalised to the numerically 
obtained result for A = 208 at y/s = 5500 GeV. The dashed curve shows the small effect of 
the term /(log(?& t /A^ CD )). 



For a thermalised system of massless bosons at an energy density e t h 
initial energy per particle can be written as 



€i, the ratio 



2.7T 



2.7^ 

U67T 2 



1/4 



E T {p s 



N( Ps 



071 Psat ^ 
\GeV) 



-0.047 



(34) 



which is indeed very close to the computed ratio E T (p Bat ) / 'N (p sa t) = £i/ n i, independent 
of A and y/s [0. Since the system looks thermal from the point of view of the average 
quantities, rapid thermalisation is plausible. This is to be contrasted with the classical 
field approach [17]] where it is found 0, |19| that the ratio Et/N is approximately 
three times larger, and therefore one might expect that thermalisation takes longer to 
achieve. On the other hand, from the analytic classical field calculation of Kovchegov 
2TJ] one infers that the e/n ratio is very close to the thermal one [2T| and, again, rapid 



thermalisation would be expected. 
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4 Final vs. initial state saturation 



Usually saturation is discussed as a small-x property of parton distribution functions. 
The above computations have been formulated referring to saturation of final state 
partons. One clearly has to understand the relation between these two approaches. 
An initial state saturation scale Q sat can naturally be defined [22| as the gluon 



transverse area density including all gluons with x > x sat = 2Q sat /y^s: 

f dyg(y,Ql t ) = -i^ocT^/GeV)^ = -QL, (35) 

■'aw. Tin, a 71 



7jR\ 7JR\ 

where the second equality was obtained by approximating xg(x, Q 2 ) = Cq{Q / GeV / x) s 
as before in Eq.(||). This equation is geometric and thus analogous to the saturation 
condition (|2l|). In the parametric weak coupling limit also this equation would contain 
various group theory factors and powers of the coupling constant. As already discussed, 
we set them equal to unity in this work. As noted below Eq. fl28|) , the powers of a s in 
the saturation condition will affect the parametric dependence of e.g. multiplicity on 

as- 
Approximating the gluon distribution as earlier and solving Q SSut from Eq. ([J5p gives 
almost the same A- and -y/i-scaling exponents as in (|22|) for p sat , only the constant is 
somewhat different and the K-f actor is absent: 

«- - ^¥-r4r^' ! ~ A " 3vf/2 ' (36) 

where Q sat and yfs are in units of GeV and Ra in fm. Note that it is essential that 
Qsat be both in the lower limit and on the right hand side of Eq.(^). The constant 
anyway is not uniquely defined, since the lower limit in (|35| ) is not unique. Thus 

Psat ~ Qsat- (37) 

In fig. (H) we plot the determination of the saturation scale using both final multiplicity 
and initial gluon distribution. From this figure we can see the small difference between 
the ^s-scaling in p sat and Q sat . 

In the initial state saturation -picture the multiplicity N of produced gluons is ex- 
pected to be proportional to N„ in Eq. (p5|), i.e. 



This relation is described in terms of the "parton liberation" constant in |p3| , and has 
been confirmed in the lattice simulations of the classical fields JL9(| . 

These results suggest that finding dynamical saturation of gluon distribution func- 
tions in a nucleus, one should also find saturation of produced gluons, the two phe- 
nomena are intimately related. 
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Figure 3: Solution of the saturation scale obtained by using the final multiplicity (solid 
lines) and Eq. (^), or the initial gluon multiplicity (dashed lines) and Eq. (35). The 
saturation scale is given by the intersection of these curves with the dotted line 'saturation' 
corresponding to PqR\- Upper two curves correspond to LHC energy while the lower two 
correspond to the full RHIC energy. Neither multiplicity contains shadowing and A = 197 



on all curves. The dashed curves could be compared with those in Fig. 2 of [22 



5 Local saturation 



In H the criterion (2~1[) was generalized to a local condition for transverse saturation 
of produced gluons in a collision with impact parameter b: 



^ = T A (s)T A (b - s)2(T pQCD (p sat ) = -p 2 sat (b, s) 



(39) 



where s is the transverse coordinate; see also [23|. Using Eq.flTE|) in ( |3"9"D one finds 
that exactly the same A and scaling exponents are obtained as from Eq. (pl~|) , and 
that the dependence on impact parameter and transverse coordinates is isolated into 
a product of nuclear density functions with a ^-dependent exponent: 



pL(b, s) ~ ^ 2 ^v^ 25/(2+?) [T A (s)T A (b - s)] 1 /^) 



and 



Mb) ~ iifVM/M / d 2 s[T A (s)T A (b-s)] 

7T J 



V(2+?) 



(40) 



(41) 



Again, 5 = 0.5 and the parameter £ as given by Eq.flTTD reproduce the behaviour of 
N(b) obtained in the numerical computation in ||. 
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With our ansatz (f|) for xg(x, Q 2 ), and neglecting the po-dependence of /(Iogpo/AQ CD ) 
in Eq. ([16]), Eq. (|39|) can also be cast into the form 

p s 2 at (b, s) ~ K^ 2 [T A ( S )T A (b _ s )]V 2 a s xg(x,pl t (b, s)) (42) 

and, consequently, 

iV(b) ~ K 1 ' 2 J d 2 s[T A (s)T A (b - s)]V^(x,^ at ( b ,s)), (43) 
where x = 2p sat (b, s)/ y's. 



Eqs. (^2|) and (^3|) permit us to comment on the relation to || , where it was postulated 
that the average (over s) saturation scale Q 2 (at fixed b) be proportional to the average 
(over s) transverse density of participating nucleons, 

Q 2 , = ^f§^,QD P ^f. (44) 

This leads to a total multiplicity, at fixed b, 

= cN pa , Tt xg(x, Q 2 S ), (45) 

with c>l, x = 2Q s /^/s. Also, the quantity N/N pait is a slightly increasing function 
of iVpart due to assumed scale evolution of the gluon structure function of the type 
xg(x,Q 2 ) ~ ln(Q s /AQCD)- The difference between and [§ can be traced down to 
two points: First, to a slightly different dependence of the saturation scale on the 
transverse coordinate originating from p par t ^ 

[T A {s)T A (h - s)] 1 / 2 . Second, to a dif- 
ferent order of averaging to obtain N(b), namely N paxt xg(x, Q 2 ) <-> / d 2 s[T A (s)T A (h — 
s)] 1 / 2 a,^(a;,p 2 at (b,s)). 



6 Discussion 

We have here shown how the A- and i/i-scaling exponents and the overall magnitude 
of various global quantities in ultrarelativistic A + A collisions, numerically computed 
in 0, can be simply related to two parameters, 5 and £. The former (Eq.(f|)) is related 
to the x~ s behaviour at small-x of the gluon distribution function at an effectively 
x-dependent saturation scale. Due to the interdependence of x and Q this is not the 
standard BFKL exponent describing small- a; behaviour at fixed scale Q. The parameter 
£ (Eq.flTTD) approximates a complicated function containing a s by a power. All the 
exponents are accurately reproduced by 5 ~ 0.5, £ = 0.44. The consequences of initial 
and final state saturation were also shown to be quantitatively similar. 
One may note the following: 
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The A-dependence of N is not that of independent hard scatterings (~ A 4 / 3 ), 
nor that of the saturation model with scaling cross section <7 p qcd ~ 1/Po (~ A) 
but, due to powerlike non-scaling of <t p qcd even slower, ~ A^ 6+2 ^^ 6+3 ^ ~ A 0,m . 
This is so even without shadowing, which further slightly reduces the exponent. A 
qualitative effect of this is that in the study of multiplicity per 0.5 times number of 
participants at some impact parameter b (which is the number to use to compare 
A+A data at various b with pp collisions) one obtains a curve decreasing very 
slowly with iV part ||]. In fact, using the simple estimate A eS = 0.5iV part (b), 
Eq . (|23|) implies that 

( £^L = \ . 0.9 • 0.8503 • K 1 ^\0.5N pai A-^ + ^^ 2S/{2+ °. (46) 
0.5iV part 3 

The decrease is thus very slow, ~ A^ 06 for f = 0.444. A more accurate 
analysis, using a local saturation condition ||, leads to a virtually constant (fa- 
independent) ratio at RHIC and slightly increasing ratio at LHC, as shown in 
Fig! 

The energy dependence of N and also of the ratio Af ch /(0.5iV part ) is the powerlike 
^■25/(2+g) ^ y^ ' 41 for 5 = 0.5, £ = 0.44. This simple power behaviour follows 
from the numerically accurate power approximation QT7j ) . This y/s dependence is 
roughly verified at RHIC for y/s = 56, 130 GeV and a new check is soon obtained 
with data at y/s = 200 GeV. At RHIC energies A+A collisions (A = 197), have 
A^ c h/(0.5A^ par t) ~ 3, clearly but not strikingly larger than the value of ~ 2 for 
p+p collisions. At LHC the increase would be from 5 for p+p to about 13 for 
A+A (Fig.£|), a really striking effect, which will directly probe the behaviour of 
the nuclear gluon densities at small values of x. 

As noted previously, the powers of a s appearing in the formulae for the mul- 
tiplicity of produced gluons, Eqs. fl28[) and fl29"[) , will be affected by additional 



powers of a s in the saturation condition (|2T|), which will appear in the weak 
coupling limit but which were replaced by constants in this study covering a 
limited energy range. However, it is interesting to note that the exponent 5 of 
the structure function appears only (at least when shadowing is neglected) in the 
\/s-scaling. The A-scaling of the multiplicity depends only on the exponent £, 
which is related to a s . Including explicitly additional powers of a s in the satura- 
tion condition (|2l|) , the £ dependence of A-scaling would change, and therefore 
the experimental measurement of A-scaling of the multiplicity would be a mea- 
surement of the actual form of the saturation criterion itself. Inclusion of, say, a 
factor 1/cKg, would make the a s dependence of multiplicity of Eq. (BBp and that 



of [13[ consistent with each other. It will also be interesting to study the relation 



to the self-screened parton cascades |2"5|. We emphasize again, however, that the 
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Figure 4: Rapidity density of charged particles near y = per 0.5 times the number of 
participants at LHC and RHIC energies computed using the local saturation criterion in |H] . 
RHIC data at y/s = 130 GeV ||] (open squares) and p+p rapidity densities (then -/V par t = 2; 
arrows) are also shown. Large (small) N par t corresponds to central b = (peripheral, b — > 
2Ra) collisions. 

purpose of this paper was to understand the scaling laws obtained numerically in 
]/J, where no explicit powers of a s were considered in the saturation condition. 
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